A Monte Carlo simulation-based decision support system for radiation oncologists in the treatment of glioblastoma multiforme.

In the present research, we have developed a model-based crisp logic function statistical classifier decision support system supplemented with treatment planning systems for radiation oncologists in the treatment of glioblastoma multiforme (GBM). This system is based on Monte Carlo radiation transport simulation and it recreates visualization of treatment environments on mathematical anthropomorphic brain (MAB) phantoms. Energy deposition within tumour tissue and normal tissues are graded by quality audit factors which ensure planned dose delivery to tumour site thereby minimising damages to healthy tissues. The proposed novel methodology predicts tumour growth response to radiation therapy from a patient-specific medicine quality audit perspective. Validation of the study was achieved by recreating thirty-eight patient-specific mathematical anthropomorphic brain phantoms of treatment environments by taking into consideration density variation and composition of brain tissues. Dose computations accomplished through water phantom, tissue-equivalent head phantoms are neither cost-effective, nor patient-specific customized and is often less accurate. The above-highlighted drawbacks can be overcome by using open-source Electron Gamma Shower (EGSnrc) software and clinical case reports for MAB phantom synthesis which would result in accurate dosimetry with due consideration to the time factors. Considerable dose deviations occur at the tumour site for environments with intraventricular glioblastoma, haematoma, abscess, trapped air and cranial flaps leading to quality factors with a lower logic value of 0. Logic value of 1 depicts higher dose deposition within healthy tissues and also leptomeninges for majority of the environments which results in radiation-induced laceration.

Dosimetric feasibility study ("proof of concept") of refractory ventricular tachycardia radioablation using proton minibeams.

PURPOSE: Preclinical data demonstrated that the use of proton minibeam radiotherapy reduces the risk of toxicity in healthy tissue. Ventricular tachycardia radioablation is an area under clinical investigation in proton beam therapy. We sought to simulate a ventricular tachycardia radioablation with proton minibeams and to demonstrate that it was possible to obtain a homogeneous coverage of an arrhythmogenic cardiac zone with this technique. MATERIAL AND METHODS: An arrhythmogenic target volume was defined on the simulation CT scan of a patient, localized in the lateral wall of the left ventricle. A dose of 25Gy was planned to be delivered by proton minibeam radiotherapy, simulated using a Monte Carlo code (TOPAS v.3.7) with a collimator of 19 0.4 mm-wide slits spaced 3mm apart. The main objective of the study was to obtain a plan ensuring at least 93% of the prescription dose in 93% of the planning target volume without exceeding 110% of the prescribed dose in the planning target volume. RESULTS: The average dose in the planning treatment volume in proton minibeam radiotherapy was 25.12Gy. The percentage of the planning target volume receiving 93% (V93%), 110% (V110%), and 95% (V95%) of the prescribed dose was 94.25%, 0%, and 92.6% respectively. The lateral penumbra was 6.6mm. The mean value of the peak-to-valley-dose ratio in the planning target volume was 1.06. The mean heart dose was 2.54Gy versus 5.95Gy with stereotactic photon beam irradiation. CONCLUSION: This proof-of-concept study shows that proton minibeam radiotherapy can achieve a homogeneous coverage of an arrhythmogenic cardiac zone, reducing the dose at the normal tissues. This technique, ensuring could theoretically reduce the risk of late pulmonary and breast fibrosis, as well as cardiac toxicity as seen in previous biological studies in proton minibeam radiotherapy.

IOeRT conventional and FLASH treatment planning system implementation exploiting fast GPU Monte Carlo: The case of breast cancer.

Partial breast irradiation for the treatment of early-stage breast cancer patients can be performed by means of Intra Operative electron Radiation Therapy (IOeRT). One of the main limitations of this technique is the absence of a treatment planning system (TPS) that could greatly help in ensuring a proper coverage of the target volume during irradiation. An IOeRT TPS has been developed using a fast Monte Carlo (MC) and an ultrasound imaging system to provide the best irradiation strategy (electron beam energy, applicator position and bevel angle) and to facilitate the optimisation of dose prescription and delivery to the target volume while maximising the organs at risk sparing. The study has been performed in silico, exploiting MC simulations of a breast cancer treatment. Ultrasound-based input has been used to compute the absorbed dose maps in different irradiation strategies and a quantitative comparison between the different options was carried out using Dose Volume Histograms. The system was capable of exploring different beam energies and applicator positions in few minutes, identifying the best strategy with an overall computation time that was found to be completely compatible with clinical implementation. The systematic uncertainty related to tissue deformation during treatment delivery with respect to imaging acquisition was taken into account. The potential and feasibility of a GPU based full MC TPS implementation of IOeRT breast cancer treatments has been demonstrated in-silico. This long awaited tool will greatly improve the treatment safety and efficacy, overcoming the limits identified within the clinical trials carried out so far.

Assessment and validation of glottic motion using cone-beam CT and real-time cine MRI.

PURPOSE: This study aimed to assess the margin for the planning target volume (PTV) using the Van Herk formula. We then validated the proposed margin by real-time magnetic resonance imaging (MRI). METHODS: An analysis of cone-beam computed tomography (CBCT) data from early glottic cancer patients was performed to evaluate organ motion. Deformed clinical target volumes (CTV) after rigid registration were acquired using the Velocity program (Varian Medical Systems, Palo Alto, CA, USA). Systematic (Σ) and random errors (σ) were evaluated. The margin for the PTV was defined as 2.5 Σ + 0.7 σ according to the Van Herk formula. To validate this margin, we accrued healthy volunteers. Sagittal real-time cine MRI was conducted using the ViewRay system (ViewRay Inc., Oakwood Village, OH, USA). Within the obtained sagittal images, the vocal cord was delineated. The movement of the vocal cord was summed up and considered as the internal target volume (ITV). We then assessed the degree of overlap between the ITV and the PTV (vocal cord plus margins) by calculating the volume overlap ratio, represented as (ITV∩PTV)/ITV. RESULTS: CBCTs of 17 early glottic patients were analyzed. Σ and σ were 0.55 and 0.57 for left-right (LR), 0.70 and 0.60 for anterior-posterior (AP), and 1.84 and 1.04 for superior-inferior (SI), respectively. The calculated margin was 1.8 mm (LR), 2.2 mm (AP), and 5.3 mm (SI). Four healthy volunteers participated for validation. A margin of 3 mm (AP) and 5 mm (SI) was applied to the vocal cord as the PTV. The average volume overlap ratio between ITV and PTV was 0.92 (range 0.85-0.99) without swallowing and 0.77 (range 0.70-0.88) with swallowing. CONCLUSION: By evaluating organ motion by using CBCT, the margin was 1.8 (LR), 2.2 (AP), and 5.3 mm (SI). The margin acquired using CBCT fitted well in real-time cine MRI. Given that swallowing during radiotherapy can result in a substantial displacement, it is crucial to consider strategies aimed at minimizing swallowing and related motion.

Assessment of bias in scoring of AI-based radiotherapy segmentation and planning studies using modified TRIPOD and PROBAST guidelines as an example.

BACKGROUND AND PURPOSE: Studies investigating the application of Artificial Intelligence (AI) in the field of radiotherapy exhibit substantial variations in terms of quality. The goal of this study was to assess the amount of transparency and bias in scoring articles with a specific focus on AI based segmentation and treatment planning, using modified PROBAST and TRIPOD checklists, in order to provide recommendations for future guideline developers and reviewers. MATERIALS AND METHODS: The TRIPOD and PROBAST checklist items were discussed and modified using a Delphi process. After consensus was reached, 2 groups of 3 co-authors scored 2 articles to evaluate usability and further optimize the adapted checklists. Finally, 10 articles were scored by all co-authors. Fleiss' kappa was calculated to assess the reliability of agreement between observers. RESULTS: Three of the 37 TRIPOD items and 5 of the 32 PROBAST items were deemed irrelevant. General terminology in the items (e.g., multivariable prediction model, predictors) was modified to align with AI-specific terms. After the first scoring round, further improvements of the items were formulated, e.g., by preventing the use of sub-questions or subjective words and adding clarifications on how to score an item. Using the final consensus list to score the 10 articles, only 2 out of the 61 items resulted in a statistically significant kappa of 0.4 or more demonstrating substantial agreement. For 41 items no statistically significant kappa was obtained indicating that the level of agreement among multiple observers is due to chance alone. CONCLUSION: Our study showed low reliability scores with the adapted TRIPOD and PROBAST checklists. Although such checklists have shown great value during development and reporting, this raises concerns about the applicability of such checklists to objectively score scientific articles for AI applications. When developing or revising guidelines, it is essential to consider their applicability to score articles without introducing bias.

Using the gamma-index analysis for inter-fractional comparison of in-beam PET images for head-and-neck treatment monitoring in proton therapy: A Monte Carlo simulation study.

GOAL: In-beam Positron Emission Tomography (PET) is a technique for in-vivo non-invasive treatment monitoring for proton therapy. To detect anatomical changes in patients with PET, various analysis methods exist, but their clinical interpretation is problematic. The goal of this work is to investigate whether the gamma-index analysis, widely used for dose comparisons, is an appropriate tool for comparing in-beam PET distributions. Focusing on a head-and-neck patient, we investigate whether the gamma-index map and the passing rate are sensitive to progressive anatomical changes. METHODS/MATERIALS: We simulated a treatment course of a proton therapy patient using FLUKA Monte Carlo simulations. Gradual emptying of the sinonasal cavity was modeled through a series of artificially modified CT scans. The in-beam PET activity distributions from three fields were evaluated, simulating a planar dual head geometry. We applied the 3D-gamma evaluation method to compare the PET images with a reference image without changes. Various tolerance criteria and parameters were tested, and results were compared to the CT-scans. RESULTS: Based on 210 MC simulations we identified appropriate parameters for the gamma-index analysis. Tolerance values of 3 mm/3% and 2 mm/2% were suited for comparison of simulated in-beam PET distributions. The gamma passing rate decreased with increasing volume change for all fields. CONCLUSION: The gamma-index analysis was found to be a useful tool for comparing simulated in-beam PET images, sensitive to sinonasal cavity emptying. Monitoring the gamma passing rate behavior over the treatment course is useful to detect anatomical changes occurring during the treatment course.

Development of porous structure for broadening Bragg-peak in scanning carbon-ion radiotherapy: Monte Carlo simulation and experimental validation.

PURPOSE: The present study aimed to develop a porous structure with plug-ins (PSP) to broaden the Bragg peak width (BPW, defined as the distance in water between the proximal and distal 80% dose) of the carbon ion beam while maintaining a sharp distal falloff width (DFW, defined as the distance along the beam axis where the dose in water reduces from 80% to 20%). METHODS: The binary voxel models of porous structure (PS) and PSP were established in the Monte Carlo code FLUKA and the corresponding physical models were manufactured by 3D printing. Both experiment and simulation were performed for evaluating the modulation capacity of PS and PSP. BPWs and DFWs derived from each integral depth dose curves were compared. Fluence homogeneity of 430 MeV/u carbon-ion beam passing through the PSP was recorded by analyzing radiochromic films at six different locations downstream the PSP in the experiment. Additionally, by changing the beam spot size and incident position on the PSP, totally 48 different carbon-ion beams were simulated and corresponding deviations of beam metrics were evaluated to test the modulating stability of PSP. RESULTS: According to the measurement data, the use of PSP resulted in an average increase of 0.63 mm in BPW and a decrease of 0.74 mm in DFW compared to PS. The 2D radiation field inhomogeneities were lower than 3 % when the beam passing through a ≥ 10 cm PMMA medium. Furthermore, employing a spot size of ≥ 6 mm ensures that beam metric deviations, including BPW, DFW, and range, remain within a deviation of 0.1 mm across various incident positions. CONCLUSION: The developed PSP demonstrated its capability to effectively broaden the BPW of carbon ion beams while maintaining a sharp DFW comparing to PS. The superior performance of PSP, indicates its potential for clinical use in the future.

Study on induced radioactivity and individual dose evaluation in Gantry room for Varian ProBeam Proton Therapy System.

Proton therapy has emerged as an advantageous modality for tumor radiotherapy due to its favorable physical and biological properties. However, this therapy generates induced radioactivity through nuclear reactions between the primary beam, secondary particles, and surrounding materials. This study focuses on systematically investigating the induced radioactivity in the gantry room during pencil beam scanning, utilizing both experimental measurements and Monte Carlo simulations. Results indicate that patients are the primary source of induced radioactivity, predominantly producing radionuclides such as 11C, 13N, and 15O. Long-term irradiation primarily generates radionuclides like 22Na, 24Na, and 54Mn etc. Additionally, this study estimates the individual doses received by medical workers in the gantry room, the irradiation dose for patient escorts, and the additional dose to patients from residual radiation. Finally, the study offers recommendations to minimize unnecessary irradiation doses to medical workers, patient escorts, and patients.

GATE Monte Carlo approach to heterogeneity dose distribution in small fields used in radiation therapy.

The Accurate dosage prediction in Radiation Therapy is challenging, prompting a need for precision beyond conventional clinical Treatment Planning Systems (TPS). Monte Carlo-based methods are sought for their superior accuracy. The aim of this study is to compare dose distributions between the ACUROS algorithm and the GATE platform in various tissue densities and field sizes, focusing on smaller fields. This study was initiated with a homogeneous validation of the TrueBeam STX system, using measurements obtained from the Centre Hospitalier Interregional Edith Cavell (CHIREC) in Brussels. The validation compared dosimetric functions (Percentage Depth Dose (PDD), Dose profile (DP) and Collimator scatter fraction (CSF)) employing the GAMMA index with a 2% / 2 mm criterion tolerance. Following this, heterogeneous studies examined dose distributions between the ACUROS algorithm and the GATE platform in various tissue densities and field sizes, with a specific focus on smaller fields. Simulations were conducted using both platforms on chest phantoms with heterogeneous slabs representing bone, lung, and heart, each housing a central tumor. The impact of electronic equilibrium on tumors for different small field sizes was evaluated. Results showed a remarkable 99% agreement between measurements and GATE calculations in the homogeneous validation of the TrueBeam STX system. However, in heterogeneous studies, ACUROS consistently overestimated lung doses by up to 8% compared to GATE simulation, especially evident with a flattening filter and smaller beam sizes at density interfaces. This highlights significant dose estimation discrepancies between ACUROS and GATE, emphasizing the need for precise calculations. The findings support exploring Monte Carlo-based methods for enhanced accuracy in Radiation Therapy treatment planning.

Aleatoric and epistemic uncertainty extraction of patient-specific deep learning-based dose predictions in LDR prostate brachytherapy.

Objective.In brachytherapy, deep learning (DL) algorithms have shown the capability of predicting 3D dose volumes. The reliability and accuracy of such methodologies remain under scrutiny for prospective clinical applications. This study aims to establish fast DL-based predictive dose algorithms for low-dose rate (LDR) prostate brachytherapy and to evaluate their uncertainty and stability.Approach.Data from 200 prostate patients, treated with125I sources, was collected. The Monte Carlo (MC) ground truth dose volumes were calculated with TOPAS considering the interseed effects and an organ-based material assignment. Two 3D convolutional neural networks, UNet and ResUNet TSE, were trained using the patient geometry and the seed positions as the input data. The dataset was randomly split into training (150), validation (25) and test (25) sets. The aleatoric (associated with the input data) and epistemic (associated with the model) uncertainties of the DL models were assessed.Main results.For the full test set, with respect to the MC reference, the predicted prostateD90metric had mean differences of -0.64% and 0.08% for the UNet and ResUNet TSE models, respectively. In voxel-by-voxel comparisons, the average global dose difference ratio in the [-1%, 1%] range included 91.0% and 93.0% of voxels for the UNet and the ResUNet TSE, respectively. One forward pass or prediction took 4 ms for a 3D dose volume of 2.56 M voxels (128 × 160 × 128). The ResUNet TSE model closely encoded the well-known physics of the problem as seen in a set of uncertainty maps. The ResUNet TSE rectum D2cchad the largest uncertainty metric of 0.0042.Significance.The proposed DL models serve as rapid dose predictors that consider the patient anatomy and interseed attenuation effects. The derived uncertainty is interpretable, highlighting areas where DL models may struggle to provide accurate estimations. The uncertainty analysis offers a comprehensive evaluation tool for dose predictor model assessment.

Multi-centre evaluation of variation in cumulative dose assessment in reirradiation scenarios.

BACKGROUND AND PURPOSE: Safe reirradiation relies on assessment of cumulative doses to organs at risk (OARs) across multiple treatments. Different clinical pathways can result in inconsistent estimates. Here, we quantified the consistency of cumulative dose to OARs across multi-centre clinical pathways. MATERIAL AND METHODS: We provided DICOM planning CT, structures and doses for two reirradiation cases: head & neck (HN) and lung. Participants followed their standard pathway to assess the cumulative physical and EQD2 doses (with provided α/ß values), and submitted DVH metrics and a description of their pathways. Participants could also submit physical dose distributions from Course 1 mapped onto the CT of Course 2 using their best available tools. To assess isolated impact of image registrations, a single observer accumulated each submitted spatially mapped physical dose for every participating centre. RESULTS: Cumulative dose assessment was performed by 24 participants. Pathways included rigid (n = 15), or deformable (n = 5) image registration-based 3D dose summation, visual inspection of isodose line contours (n = 1), or summation of dose metrics extracted from each course (n = 3). Largest variations were observed in near-maximum cumulative doses (25.4 - 41.8 Gy for HN, 2.4 - 33.8 Gy for lung OARs), with lower variations in volume/dose metrics to large organs. A standardised process involving spatial mapping of the first course dose to the second course CT followed by summation improved consistency for most near-maximum dose metrics in both cases. CONCLUSION: Large variations highlight the uncertainty in reporting cumulative doses in reirradiation scenarios, with implications for outcome analysis and understanding of published doses. Using a standardised workflow potentially including spatially mapped doses improves consistency in determination of accumulated dose in reirradiation scenarios.

Dosimetric evaluation in Helical TomoTherapy for lung SBRT using Monte Carlo-based independent dose verification software.

PURPOSE: To elucidate the dosimetric errors caused by a model-based algorithm in lung stereotactic body radiation therapy (SBRT) with Helical TomoTherapy (HT) using Monte Carlo (MC)-based dose verification software. METHODS: For 38 plans of lung SBRT, the dose calculation accuracy of a treatment planning system (TPS) of HT was compared with the results of DoseCHECK, the commercial MC-based independent verification software. The following indices were extracted to evaluate the correlation of dosimetric errors: (1) target volume, (2) average computed tomography (CT) value of the planning target volume (PTV) margin, and (3) average CT value of surrounding 2-mm area of the PTV (PTV ring). Receiver operating characteristic (ROC) analyses determined the threshold for 5% of differences in PTV D95%. Then, the 38 plans were classified into two groups using the cutoff values of ROC analysis for these three indices. Dosimetric differences between groups were statistically compared using the Mann-Whitney U test. RESULTS: TPS of HT overestimated by more than 5% in the PTV D95% in 16 of 38 plans. The PTV ring showed the strongest correlation with dosimetric differences. The cutoff value for the target volume, the PTV margin, and the PTV ring was 14.7 cc, -754 HU, and -708 HU, respectively. The area under the curve (AUC) for the target volume, the PTV margin, and the PTV ring were 0.835, 0.878, and 0.932, respectively. Dosimetric errors more than 5% were observed when the PTV volume was less than 15 cc or when the CT value around the target was less than -700 HU. CONCLUSION: The TPS of HT might overestimate the PTV dose by more than 5% if any the three indices in this study were below threshold. Therefore, independent verification with an MC-based algorithm should be strongly recommended for lung SBRT in HT.

Very high-energy electron therapy as light-particle alternative to transmission proton FLASH therapy - An evaluation of dosimetric performances.

PURPOSE: Clinical translation of FLASH-radiotherapy (RT) to deep-seated tumours is still a technological challenge. One proposed solution consists of using ultra-high dose rate transmission proton (TP) beams of about 200-250 MeV to irradiate the tumour with the flat entrance of the proton depth-dose profile. This work evaluates the dosimetric performance of very high-energy electron (VHEE)-based RT (50-250 MeV) as a potential alternative to TP-based RT for the clinical transfer of the FLASH effect. METHODS: Basic physics characteristics of VHEE and TP beams were compared utilizing Monte Carlo simulations in water. A VHEE-enabled research treatment planning system was used to evaluate the plan quality achievable with VHEE beams of different energies, compared to 250 MeV TP beams for a glioblastoma, an oesophagus, and a prostate cancer case. RESULTS: Like TP, VHEE above 100 MeV can treat targets with roughly flat (within ± 20 %) depth-dose distributions. The achievable dosimetric target conformity and adjacent organs-at-risk (OAR) sparing is consequently driven for both modalities by their lateral beam penumbrae. Electron beams of 400[500] MeV match the penumbra of 200[250] MeV TP beams and penumbra is increased for lower electron energies. For the investigated patient cases, VHEE plans with energies of 150 MeV and above achieved a dosimetric plan quality comparable to that of 250 MeV TP plans. For the glioblastoma and the oesophagus case, although having a decreased conformity, even 100 MeV VHEE plans provided a similar target coverage and OAR sparing compared to TP. CONCLUSIONS: VHEE-based FLASH-RT using sufficiently high beam energies may provide a lighter-particle alternative to TP-based FLASH-RT with comparable dosimetric plan quality.

A time- and space-saving Monte Carlo simulation method using post-collimation generative adversarial network for dose calculation of an O-ring gantry Linac.

PURPOSE: This study explores the feasibility of employing Generative Adversarial Networks (GANs) to model the RefleXion X1 Linac. The aim is to investigate the accuracy of dose simulation and assess the potential computational benefits. METHODS: The X1 Linac is a new radiotherapy machine with a binary multi-leaf collimation (MLC) system, facilitating innovative biology-guided radiotherapy. A total of 34 GAN generators, each representing a desired MLC aperture, were developed. Each generator was trained using a phase space file generated underneath the corresponding aperture, enabling the generation of particles and serving as a beam source for Monte Carlo simulation. Dose distributions in water were simulated for each aperture using both the GAN and phase space sources. The agreement between dose distributions was evaluated. The computational time reduction from bypassing the collimation simulation and storage space savings were estimated. RESULTS: The percentage depth dose at 10 cm, penumbra, and full-width half maximum of the GAN simulation agree with the phase space simulation, with differences of 0.4 % ± 0.2 %, 0.32 ± 0.66 mm, and 0.26 ± 0.44 mm, respectively. The gamma passing rate (1 %/1mm) for the planar dose exceeded 90 % for all apertures. The estimated time-saving for simulating an plan using 5766 beamlets was 530 CPU hours. The storage usage was reduced by a factor of 102. CONCLUSION: The utilization of the GAN in simulating the X1 Linac demonstrated remarkable accuracy and efficiency. The reductions in both computational time and storage requirements make this approach highly valuable for future dosimetry studies and beam modeling.

Neutron flux evaluation algorithm with a combination of Monte Carlo and removal-diffusion calculation methods for boron neutron capture therapy.

BACKGROUND: In Japan, the clinical treatment of boron neutron capture therapy (BNCT) has been applied to unresectable, locally advanced, and recurrent head and neck carcinomas using an accelerator-based neutron source since June of 2020. Considering the increase in the number of patients receiving BNCT, efficiency of the treatment planning procedure is becoming increasingly important. Therefore, novel and rapid dose calculation algorithms must be developed. We developed a novel algorithm for calculating neutron flux, which comprises of a combination of a Monte Carlo (MC) method and a method based on the removal-diffusion (RD) theory (RD calculation method) for the purpose of dose calculation of BNCT. PURPOSE: We present the details of our novel algorithm and the verification results of the calculation accuracy based on the MC calculation result. METHODS: In this study, the "MC-RD" calculation method was developed, wherein the RD calculation method was used to calculate the thermalization process of neutrons and the MC method was used to calculate the moderation process. The RD parameters were determined by MC calculations in advance. The MC-RD calculation accuracy was verified by comparing the results of the MC-RD and MC calculations with respect to the neutron flux distributions in each of the cubic and head phantoms filled with water. RESULTS: Comparing the MC-RD calculation results with those of MC calculations, it was found that the MC-RD calculation accurately reproduced the thermal neutron flux distribution inside the phantom, with the exception of the region near the surface of the phantom. CONCLUSIONS: The MC-RD calculation method is useful for the evaluation of the neutron flux distribution for the purpose of BNCT dose calculation, except for the region near the surface.

Systematic Implementation of Effective Quality Assurance Processes for the Assessment of Radiation Target Volumes in Head and Neck Cancer.

PURPOSE: Significant heterogeneity exists in clinical quality assurance (QA) practices within radiation oncology departments, with most chart rounds lacking prospective peer-reviewed contour evaluation. This has the potential to significantly affect patient outcomes, particularly for head and neck cancers (HNC) given the large variance in target volume delineation. With this understanding, we incorporated a prospective systematic peer contour-review process into our workflow for all patients with HNC. This study aims to assess the effectiveness of implementing prospective peer review into practice for our National Cancer Institute Designated Cancer Center and to report factors associated with contour modifications. METHODS AND MATERIALS: Starting in November 2020, our department adopted a systematic QA process with real-time metrics, in which contours for all patients with HNC treated with radiation therapy were prospectively peer reviewed and graded. Contours were graded with green (unnecessary), yellow (minor), or red (major) colors based on the degree of peer-recommended modifications. Contours from November 2020 through September 2021 were included for analysis. RESULTS: Three hundred sixty contours were included. Contour grades were made up of 89.7% green, 8.9% yellow, and 1.4% red grades. Physicians with >12 months of clinical experience were less likely to have contour changes requested than those with <12 months (8.3% vs 40.9%; P < .001). Contour grades were significantly associated with physician case load, with physicians presenting more than the median number of 50 cases having significantly less modifications requested than those presenting <50 (6.7% vs 13.3%; P = .013). Physicians working with a resident or fellow were less likely to have contour changes requested than those without a trainee (5.2% vs 12.6%; P = .039). Frequency of major modification requests significantly decreased over time after adoption of prospective peer contour review, with no red grades occurring >6 months after adoption. CONCLUSIONS: This study highlights the importance of prospective peer contour-review implementation into systematic clinical QA processes for HNC. Physician experience proved to be the highest predictor of approved contours. A growth curve was demonstrated, with major modifications declining after prospective contour review implementation. Even within a high-volume academic practice with subspecialist attendings, >10% of patients had contour changes made as a direct result of prospective peer review.

Plan quality assessment of modern radiosurgery technologies in the treatment of multiple brain metastases.

Stereotactic radiosurgery (SRS) of multiple brain metastases has evolved over the last 40 years allowing centres to treat an increasing number of brain metastases in a single treatment fraction. HyperArcTMplanning optimisation technique is one such development that streamlines the treatment of multiple metastases with a single isocentre. Several studies have investigated the plan quality of HyperArc compared to CyberKnife or Gamma Knife, however there are limited number of studies that include all three modalities. It is the aim of this study to provide an assessment of plan quality between the three SRS platforms across ten patients with multiple brain metastases ranging from three to eight metastases per patient. Strict planning workflows were established to avoid bias towards any particular treatment platform. Plan quality was assessed through dose to organs at risk, Paddick conformity index (PCI), gradient index (GI), global efficiency index (Gη) and dose to normal brain tissue. Results from this study found mean PCI observed across Gamma Knife plans was significantly lower than HyperArc and CyberKnife. HyperArc plans observed significantly shorter beam-on times which were 10 to 20 times faster than CyberKnife and Gamma Knife plans. Gamma Knife and CyberKnife were found to produce plans with significantly superior GI, global efficiency index and the volume of healthy brain receiving greater than 12 Gy (V12Gy) when compared to HyperArc plans. Lesion volume was seen to influence the relative difference in dose metrics between systems. The study revealed that all three treatment modalities produced high quality plans for the SRS treatment of multiple brain metastases, each with respective benefits and limitations.

Development and optimisation of grid inserts for a preclinical radiotherapy system and corresponding Monte Carlo beam simulations.

Objective. To develop a physical grid collimator compatible with the X-RAD preclinical radiotherapy system and create a corresponding Monte Carlo (MC) model.Approach. This work presents a methodology for the fabrication of a grid collimator designed for utilisation on the X-RAD preclinical radiotherapy system. Additionally, a MC simulation of the grid is developed, which is compatible with the X-RAD treatment planning system. The grid was manufactured by casting a low melting point alloy, cerrobend, into a silicone mould. The silicone was moulded around a 3D-printed replica of the grid, enabling the production of diverging holes with precise radii and spacing. A MC simulation was conducted on an equivalent 3D grid model and validated using 11 layers of GAFChromic EBT-3 film interspersed in a 3D-printed water-equivalent phantom. A 3D dose distribution was constructed from the film layers, enabling a direct comparison with the MC Simulation.Main results. The film and the MC dose distribution demonstrated a gamma passing rate of 99% for a 1%, 0.5 mm criteria with a 10% threshold applied. The peak-to-valley dose ratio and output factor at the surface were determined to be 20.4 and 0.79, respectively.Significance. The pairing of the grid collimator with a MC simulation can significantly enhance the practicality of grid therapy on the X-RAD. This combination enables further exploration of the biological implications of grid therapy, supported by a knowledge of the complex dose distributions. Moreover, this methodology can be adapted for use in other systems and scenarios.

Clinical assessment of deep learning-based uncertainty maps in lung cancer segmentation.

Objective. Prior to radiation therapy planning, accurate delineation of gross tumour volume (GTVs) and organs at risk (OARs) is crucial. In the current clinical practice, tumour delineation is performed manually by radiation oncologists, which is time-consuming and prone to large inter-observer variability. With the advent of deep learning (DL) models, automated contouring has become possible, speeding up procedures and assisting clinicians. However, these tools are currently used in the clinic mostly for contouring OARs, since these systems are not reliable yet for contouring GTVs. To improve the reliability of these systems, researchers have started exploring the topic of probabilistic neural networks. However, there is still limited knowledge of the practical implementation of such networks in real clinical settings.Approach. In this work, we developed a 3D probabilistic system that generates DL-based uncertainty maps for lung cancer CT segmentations. We employed the Monte Carlo (MC) dropout technique to generate probabilistic and uncertainty maps, while the model calibration was evaluated by using reliability diagrams. A clinical validation was conducted in collaboration with a radiation oncologist to qualitatively assess the value of the uncertainty estimates. We also proposed two novel metrics, namely mean uncertainty (MU) and relative uncertainty volume (RUV), as potential indicators for clinicians to assess the need for independent visual checks of the DL-based segmentation. Main results. Our study showed that uncertainty mapping effectively identified cases of under or over-contouring. Although the overconfidence of the model, a strong correlation was observed between the clinical opinion and MU metric. Moreover, both MU and RUV revealed high AUC values in discretising between low and high uncertainty cases.Significance. Our study is one of the first attempts to clinically validate uncertainty estimates in DL-based contouring. The two proposed metrics exhibited promising potential as indicators for clinicians to independently assess the quality of tumour delineation.

Pencil beam kernel-based dose calculations on CT data for a mixed neutron-gamma fission field applying tissue correction factors.

Objective.For fast neutron therapy with mixed neutron and gamma radiation at the fission neutron therapy facility MEDAPP at the research reactor FRM II in Garching, no clinical dose calculation software was available in the past. Here, we present a customized solution for research purposes to overcome this lack of three-dimensional dose calculation.Approach.The applied dose calculation method is based on two sets of decomposed pencil beam kernels for neutron and gamma radiation. The decomposition was performed using measured output factors and simulated depth dose curves and beam profiles in water as reference medium. While measurements were performed by applying the two-chamber dosimetry method, simulated data was generated using the Monte Carlo code MCNP. For the calculation of neutron dose deposition on CT data, tissue-specific correction factors were generated for soft tissue, bone, and lung tissue for the MEDAPP neutron spectrum. The pencil beam calculations were evaluated with reference to Monte Carlo calculations regarding accuracy and time efficiency.Main results.In water, dose distributions calculated using the pencil beam approach reproduced the input from Monte Carlo simulations. For heterogeneous media, an assessment of the tissue-specific correction factors with reference to Monte Carlo simulations for different tissue configurations showed promising results. Especially for scenarios where no lung tissue is present, the dose calculation could be highly improved by the applied correction method.Significance.With the presented approach, time-efficient dose calculations on CT data and treatment plan evaluations for research purposes are now available for MEDAPP.